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Abstract 

The spontaneous creation and persistence of ground-state coherence in an ensemble of intracav- 
ity Rb atoms has been observed as a quantum beat. Our system realizes a quantum eraser, where 
the detection of a first photon prepares a superposition of ground-state Zeeman sublevels, while de- 
tection of a second erases the stored information. Beats appear in the time-delayed photon-photon 
coincidence rate (intensity correlation function). We study the beats theoretically and experimen- 
tally as a function of system parameters, and find them remarkably robust against perturbations 
such as spontaneous emission. Although beats arise most simply through single-atom-mediated 
quantum interference, scattering pathways involving pairs of atoms interfere also in our intracavity 
experiment. We present a detailed model which identifies all sources of interference and accounts 
for experimental realities such as imperfect pre-pumping of the atomic beam, cavity birefringence, 
and the transit of atoms across the cavity mode. 

PACS numbers: 42.50.Md,42.50.Pq,37.30.+i 
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I. INTRODUCTION 



The interference of scattering amplitudes in quantum mechanics arises from the indistin- 
guishability of alternative scattering paths, as in the interference of the paths of a photon 
passing through the celebrated double-slits of Young. For example, a modern variation 
on the Young experiment jlj] shows spatial fringes in the intensity of light scattered from 
two trapped ions; but only so long as the scattering cannot be traced to one ion or the 
other [2j. Even in the latter situation, other measurements might be found which cannot 
distinguish between the paths, hence recovering the interference. Typically they involve 
higher-order moments of the field jsj, or post-selection, i.e., partitioning of scattering events 
into subensembles jj]. Such a measurement strategy is termed a "quantum eraser", since it 
recovers interference by "erasing" the information that identifies the path. 

A time-domain analogue of double-slit interference can occur inside multi-level atoms, 
where photon emission via parallel transitions can result in a modulation of the emission 
intensity at the frequency of a level splitting — the phenomenon of "quantum beats" (see, 
e.g., 0, 0]). A distinction has traditionally been drawn between "Type-I" or "V" atomic 
systems, where decay of a superposition of upper levels yields beats at the transition differ- 
ence frequency; and the inverted "Tvpe-II" or "A" systems, where decay to a superposition 
of lower levels does not yield beats [3, |8] . The typical argument for the latter outcome is 
that a measurement of the ground-state population could always, in principle, determine in 
which of the two available states the electron landed; as there is no sum over alternative 
paths to one and the same final state, there are no "ground-state quantum beats" 0, Q . 

Nevertheless, as with spatial fringes, a quantum eraser-type strategy can recover time- 
domain interference in the ground state. Zajonc Q, [ll| proposed one such implementation 
in a "Type-II" atomic system, basing his proposal on two-photon scattering. In this case 
the second scattered photon erases the path information written by the first — amplitudes 
for the scattering of two photons in sequence interfere. 

We recently published experimental results showing quantum beats in spontaneous emis- 
sion at the frequency of the ground-state Zeeman splitting in Rb, i.e. ground-state quantum 



beats seen in spontaneous emission 



121 ] . Oscillations appear in the second-order intensity 



autocorrelation function only, not in the average intensity, as follows from the indistinguisha- 
bility requirement above. The presence of a similar oscillation hidden within the noise of 
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spontaneous emission was demonstrated in 1955 by Forrester et al. [13J. They mixed two 
incoherent light sources — a Zeeman doublet — on a photocathode, and used a resonant mi- 
crowave cavity to enhance the beat signal extracted from the photocurrent noise. The 
interference in this case is classical, though the oscillation is recovered from noise through 
intensity correlation. More than 50 years later, using coherent excitation and single-photon 
detectors, we have realized a time-resolved measurement of the ground-state quantum beat 
recovered from spontaneous emission noise. 

The oscillation in our system arises as a complicated mixture of quantum eraser-type 
interferences within single atoms, pair-wise interference between emission from different 
atoms, and a homo dyne contribution due to the superposition of a weak coherent back- 
ground (similar to 14|) generated by birefringence in the cavity mirrors. By coupling spon- 
taneous emission into an optical cavity at moderate dipole coupling stren gth, we overcome 
the signal-to-noise limitations set by a small coherence area in free space 13j, and enforce 
indistinguishability among different atoms emitting into a common spatial mode. Moreover, 
we show below that the complicated level structure of 85 Rb actually aids in the survival 
of ground-state coherence, counter to the conventional strategy of protecting coherence by 
limiting the state space through which population can diffuse. 

We distinguish here between ground-state coherences imposed by an external drive and 
those arising spontaneously, selected through the detection process, as in our experiment. 
In the former case, an external magnetic or optical drive couples two ground states directly, 
with the resulting coherence read out optically in forward scattering (see [jjj for many 
examples.) In the latter, levels couple only through the vacuum, with no external drive 
to enforce coherence. The fact that spontaneous decay can generate coherence is evident 
rom the observation of quantum beats at the intermediate level splitting in cascade decay 
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171 ]; that the same process occurs in transitions to ground or meta-stable states is not 



therefore surprising. Schubert et al. [18] measured such a coherence in the bichromatic 
cross-correlation of fluorescence from a single ion, where detection of a first photon left the 
ion in a superposition of meta-stable states. When considering isotropic emission, however, 
spontaneously created coherences tend to vanish on the average, and for this reason are 
often left out of density matrix calculations 19j. The recent interest follows a 1992 paper 



20j in which specific measurab 



been explored theoretically 21 



e consequences were claimed; various arrangements have 



22j. A publication as recent as 2005 23J claims evidence 
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of the first serious experimental consequences of spontaneously-generated coherence, this in 
a quantum dot system. Other experiments are surprisingly few. We direct the interested 
reader to Ref. 24( for an overview. 

In this article we expand upon the results presented in [l^ . In particular, we seek to 
explain the origin of the various individual components of the beat signal, show which 
experimental conditions are necessary for the robust survival and detection of beats, and 
explore their sensitivity to various experimental controls. The paper is organized as follows. 
Section [III introduces the theoretical model, starting with a single atom fixed in space and 
moving to a full atomic beam with realistic fluctuations. Section IIHI presents the details of 
our experimental method, from the atomic source and optics to the detection apparatus. 
Section [111 Bl summarizes the evolution of the beat signal as we explore parameter space and 
compares experiment with theory. The paper concludes in Section IIVI 



II. THEORETICAL MODEL 

We consider first the idealized system of one fixed 85 Rb atom, then turn to a realistic 
atomic ensemble, as realized in our experiment with a cold atomic beam. The atom has 
Zeeman structure in its ground and excited states [Fig. QJa)] and interacts through the 
D2 line with degenerate, orthogonally polarized cavity modes, designated H (horizontal) 
and V (vertical); a weak magnetic field sets the quantization axis in the vertical direction, 
and mode V is weakly and continuously driven [Fig. QJb)]. The atom is prepared in state 
|<7o)> from which it is excited to state |ep) by the V mode [Fig. 12(a)]- It may return to 
the ground state by emitting a it, <j + or o~ photon, or any linear combination conserving 
angular momentum. In the assigned geometry, only <j + or a~ light couples to the H mode, 
with the helicity undetermined. We assume that the probability of reabsorption of an 
emitted photon is negligible; it escapes the cavity and its detection places the atom in the 
superposition [Fig. 12(b)]: 

|^o> = (l<M> + l0+i>)/V2. (1) 

The atom is now in the ground state with its angular momentum perpendicular to the 
magnetic field, and it performs Larmor precession. With subsequent reexcitation by the V 
mode, the state 

|^(t)) = (e-^|e_ 1 ) + e^|e +1 ))/v / 2 (2) 
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FIG. 1: (Color online) Experimental setup: (a) 7r-excitation of the F = 3 to F' = 4 transition in 
Rb showing scattering of a first (red) and second (blue) photon into the H mode; (b) schematic 
of the apparatus, HWP: half-wave plate, PBS: polarizing beam-splitter, BS: beam-splitter, APD: 
avalanche photo diode. 

is reached, with phase ±0(t) gained through precession in the ground state [Fig. [2]^c)]. From 
here the atom can decay back to \go) by emitting a second if -mode photon [Fig. EJ^d)]. The 
probability for this emission depends on the phase difference, 2<p(t), between the two parts 
of the superposition. It oscillates and thus gives rise to beats in the rate of detection of a 
second H photon subsequent to the detection of a first. 

As depicted in Fig. [2], there are two paths for scattering a pair of photons into the H 
mode: \g Q ) -> |e ) -> \g+i) ->■ |e +1 ) ->• \g ) and \g Q ) ->■ |e ) ->■ \g_ t ) -» |e_i> ->■ \g ). The 
phase gained from the ground-state Zeeman shift (Larmor precession) is different along the 
two paths, which interfere to produce oscillations in the rate of delayed coincidences — i.e., in 
the correlation function g^(r). Note that after the first photon is detected, "which path" 
information is available, since \g+i) and \g-i) are distinguishable in principle, and their 
orthogonality precludes observation of interference effects in the average intensity arising 
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(a) (b) 




FIG. 2: Simplified model of the two-photon quantum eraser process, (a) 7r-excitation from go to 
eo; (b) spontaneous decay through a transitions to a superposition of g-\ and g + \; (c) 7r-excitation 
to a superposition of e_i and e+i; (d) spontaneous decay through cr transitions back to go- 

from cross-terms when taking the expectation value over Eq. [TJ This information is largely 
erased by the second photon detection, where the amplitude for returning to the common 
final state \go) allows survival of cross-terms in the two-photon intensity expectation value, 
corresponding to interference between the two scattering paths. We note that, as in Ref. [l^], 
the spectrum of scattered light still exhibits a doublet separated by the beat frequency, but 
the random phase relation between the two fields destroys any first-order coherence in the 



6 



average intensity. 

In the following sections we study the dependence of the quantum beats on the different 
features and parameters of our experimental system: magnetic field strength, number of 
atoms, initial state preparation, atomic beam fluctuations, cavity birefringence, and mixing 
of the undriven H field with drive light (V field) outside the cavity before detection. We first 
present a detailed theory for the case of one fixed atom before we move on to the treatment 
of many atoms and atomic motion. 



A. One fixed atom 



The relevant 16- level structure of the atom is depicted in Fig. [TJ Bold black arrows in- 
dicate its interaction with the driven V mode of the optical cavity, and red and blue wavy 
lines with the undriven H mode. We calculate the second-order correlation function of the 



H mode. We use the quantum trajectory formalism [25j|, which provides insight into the 
physical processes involved and facilitates efficient numerical calculations, something of im- 
portance when many atoms are considered. Working in a frame rotating at the frequency 
of the drive, the non-Hermitian Hamiltonian governing coherent evolution between spon- 
taneous emission events (modes other than H and V) or photon loss through the cavity 
mirrors is: 



with free Hamiltonian 



interaction and drive 



H S = H + Hj + H D + H L , (3) 



H = h5 v ala v + Ti5ha} h a,h 

3 4 

+ H h&Bi\9%)(Si\ + h8 ei \ei)(ei\ , (4a) 

i=— 3 i=— 4 



Hj = hg [a+E, + a\ (E ff+ + E a _ )] + h.c. , (4b) 
H D = ih£(al — a v ) , (4c) 



and non-Hermitian loss term 



Hl = —ihn(ala v + a\ah) 

+ + Et_Ea_) , (4d) 
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where a v and ah annihilate photons in the V and H modes, respectively, detunings from the 
drive 5 V and 5h, £ is the drive amplitude for the V mode, g is the dipole coupling constant, 
2k is the photon loss rate from each cavity mode, 7 is the spontaneous emission rate, and 
H n , £<j + , and E CT _ are dipole lowering operators for n, cr + , and cr_ transitions; atomic energy 
shifts in the magnetic field and relative to the drive are HS ei , i = —4, —3, ... ,4 (excited 
states) and H5 gi , i = —3, —2, . . . , 3 (ground states). Explicit forms for the dipole operators 
depend on Clebsh-Gordon coefficients and are relegated to the Appendix. 

The evolution of the system under the non-Hermitian Hamiltonian H$ is calculated nu- 
merically, with photon-number truncation, using the fourth-order Runge-Kutta method, 
including step-size correction. At regular intervals of length At the atom and the cavity 
modes are checked for quantum jumps — spontaneous emission or photon leakage. If a scat- 
tering event occurs, the system state is collapsed accordingly — jump operators E^, E CT+ , 
E(j_, a v , or ah — before continuing the coherent evolution. After a time t, we assume a 
photon is scattered through the H mode. The system state is collapsed (jump operator 
ah) and its conditional evolution, with the unnormalized correlation function, G^(t, t + r), 
given by the i7-mode-photon-number-weighted mean of the i7-mode photon number ex- 
pectation in the conditional states. Normalization by the photon number averages yields 
g W(t, t + r)= G<?){t, t + r)/((ala h )(t))((aia h )(t + r)>. 

We consider now a weak drive, such that (a\a v ) <C 1, and t long enough for the system 
to reach a quasi-steady-state (overlooking the very slow process of optical pumping). If the 
initial state is chosen to be any one of the ground states \gi), i = — 2, — 1, 0, +1, +2, the 
correlation function, g^(r), shows beats, with beat frequency twice the Larmor precession 
frequency of \g±i); the frequency does not depend on the chosen initial state, or the number 
of atoms when many are included (Sec. Ill Bp . 

The visibility of the beat depends, however, on a number of issues. In the ideal case, the 
atom is in state \go) when the first if-photon is detected; it is projected to |-?/>o), and the 
second photon is detected as it returns to \go) via \ipb(t)) (Fig. [2]). We note first that, for 
the level scheme of Fig. (TJ the return is not to \go) but (see the Appendix for the explicit 
form of the lowering operator, Y} a+ + E^ , for if-mode scattering) 

2cos[ < Kt)]^ > + \l^(e^%-2) + e^\g + ,)). (5) 
This brings a reduction of the visibility to a little more than 75%. Beyond this, several 
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Initial state 


9o 9+1 9+2 9+3 


Visibility 


0.75 0.5 0.15 0.03 



TABLE I: Quantum beat visibility for one fixed atom and different initial ground states. The 
parameters are: g/'j = 0.25, re/7 = 0.5, £/j = 0.025, and a Larmor frequency of 2.27/3. The time 
at which the first H-mode photon is detected is t = 257 -1 . Results for g_i, g_2, g-3 follow by 
symmetry. Note that a visibility of zero (no quantum beat) is predicted for g ± 3 in the absence of 
optical pumping prior to t. 

trajectories deviate from the ideal and further reduce the visibility. We divide them into 
two groups: (i) those that deviate prior to the detection of the first photon, and (ii) those 
that deviate between the detection of the first photon and the second. 

Consider first group (i). If the atom is initially in state \gi), % 7^ 0, then when the first H- 
photon is detected, it is projected into an unequal superposition of \gi±i), with the different 
weights given by Clebsch-Gordan coefficients; for example, if the initial state is \g+i), after 
detection of a scattered photon, the atom is projected into \ip) = (vT0|g ) + V / 3|5'2))/"\/l3- 
Table [T] lists the realized visibilities for all potential initial states. A maximum visibility 
of 0.75 is achieved when the initial state is \go). The visibility decays rather quickly away 
from this maximum. Experimentally, we optically pump the atoms before they enter the 
cavity; however, the efficiency of the pumping is not perfect, and the result is generally a 
distribution over ground states, peaked around \go). Some reduction of the beat visibility 
must follow from the imperfect optical pumping prior to an atom entering the cavity. 

Other processes can redistribute population amongst the atomic levels and contribute 
to the loss of visibility. The atom can spontaneously decay from an excited state |e$) to 
|<7i+i) or |<7i_i) with the emission of a photon to the side rather than into the H mode 
of the cavity. At higher values of the drive, this process might be repeated many times, 
redistributing population before a first if-photon is detected. The distribution reached 
through such optical pumping by the drive depends on the Clebsch-Gordan coefficients, 
the drive strength, and the time when the first photon is detected. Table [TT] displays the 
distribution reached in the long-time limit. It shows that, even in the limit, states with 
i = ±3 have a very small probability to be populated, and those with i = ±2 are populated 
at a level of only about 10%. 
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The distribution over atomic ground states prior to the detection of a first if-photon 
strongly effects the state the atom is projected into on average. If the time of the first 
detection is close to zero the atom is projected into [Eq. ([T]. As this time increases, 
superpositions of, first, \go) and \g+ 2 ) (or |^r_ 2 )), and then \g±i) and \g±s) appear. Numeri- 
cally we have checked that the state immediately after the first if -photon detection may be 
written approximately as 

Pc(t) = Po\lpo)(lpo\+Pl(\^l){H + \^-l}(^-l\) 

+p 2 Ap c , (6) 

with 

|V±i) = (Vl0\g } + V3\g ±2 })/Vl3, (7) 

where pi, i = 0, 1, is the probability distribution over \ip ) and \ip±i), and p 2 = 1 — Po — 2pi 
is the probability that the first //-photon is scattered out of one of the ground states with 
\i\ = 2,3; Ap c is the state reached from such scattering events. As can be seen from 
Fig. [31 as the wait for the first if-photon detection becomes longer, the probability to 
realize \ip s ) decreases while that to realize l^ii) increases. For the parameters considered, 
the sum p + 2pi is close to 0.97, which tells us that \ip s ) and l^ii) cover all relevant 
prepared superpositions. An increase in drive strength changes only the time evolution of 
the probabilities, not their stationary values. 

We move now to trajectories that deviate from the ideal after the first //-photon is 
detected [group (ii)]. We assume the prepared superposition is |^o). Interaction of the atom 
with the driven V cavity mode moves population to the superposition \ip' Q (t)) [Eq. (J2J)]. It is 
then possible that spontaneous emission (to the side) moves \i[>b(t)) to a superposition of the 
ground states \go) and \g+ 2 ) (or \g- 2 ))- As the Clebsch-Gordan coefficient connecting |e T i) 
to \g ) differs from that connecting |e±i) to ^±2)), such an event yields unequal weights 



State 

Probability 



9o 9+1 5+2 5+3 e e +1 e +2 e +3 
0.23 0.15 0.04 0.003 0.16 0.09 0.02 0.001 



TABLE II: Energy level occupation probabilities for one fixed atom. The parameters are: 5/7 = 
0.25, re/7 = 0.5, 5/7 = 2, and a Larmor frequency of 2.27/3. The number of V-mode photons 
inside the cavity is (ala v ) 15. Results for g^ 2 , g_3 and e_i, e_2, e_3 follow by symmetry. 
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Yt 

FIG. 3: Probability for the preparation of superposition state \ip s ) (solid line) and either of the 
superposition states \ip±i) (dashed line), for one fixed atom, as a function of the time of the first 
-£/-photon detection. The parameters are: 5/7 = 0.25, re/7 = 0.5, £/7 = 0.3, and a Larmor 
frequency of 2.27/3. 

in the ground-state superposition. Nevertheless, this aside, the described process recovers 
the initial setup — a ground-state superposition — but in a manifold of states shifted to the 
right or left. Continuing then with the standard story, the 7r-polarized drive transfers the 
superposition of \g ) and ^±2) to a superposition of |e ) and |e± 2 ) in the excited state. A 
second if -photon can then be emitted via the cavity, projecting the atom into a superposition 
°f |#±i), ^±3) an d |<7=fi) [compare Eq. (|5])]. Apart for the changed weight factors, the 
quantum eraser process still takes place, only within a different manifold of atomic states. 
Of course, the unequal weights yield a beat with diminished visibility. 

Such spontaneous emission events can happen several times in the interval separating the 
ii-photon detections; nevertheless, so long as there is no a + (o~_) emission at a time when 
|e_3) (|e + 3)) is part of the superposition in the excited state, the superposition, with modified 
weights, survives. For an atom with six levels, as in Fig. ((2]), a single cr+ or er_ spontaneous 
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emission will destroy the prepared superposition, reducing the observed visibility far more 
than in the 16-level case. Our use of the 16-level configuration produces a particularly robust 
(against spontaneous emission) quantum beat. 

Figure H] shows how the spontaneously created coherence moves between different ground 
states after the first .£/-photon is detected. We quantify the coherence by off-diagonal matrix 
elements \(gi+i\p c (t + r)|<7j_i)|, i = and 1, which fall between a maximum of 0.5 (equally 
weighted superposition) and zero. The solid line follows | (g + i\p(t-\-r) | | as a function of r. 
It begins at r = from approximately 0.46, which shows that immediately after the photon 
is detected the atom is to a good approximation in the superposition \^o); po ~ 1 in Eq. (|6]). 
As the time to the second photon detection progresses, p decreases as p±i grows and part of 
the coherence is transferred to a superposition of \g ) and \g±2)- Eventually phase diffusion 
(decoherence) which accompanies repeated cycles of excitation and spontaneous emission 



sets in, causing both displayed coherences to decay to zero [26J . Coherences between ground 
states other than those shown in the figure are negligible. 

Before turning to the many atom case relevant to our experiment, one final effect might 
usefully be introduce at the one-atom level. The mirrors in the experiment show a small 
birefringence and mix a little of the ^-polarized light with the if-polarized cavity mode. We 
attempt to null this mixing with a half- wave plate placed in the cavity output [Fig. Q^b)]; 
alternatively, in some measurements we deliberately enhance it (see Fig. [TT]) . The mixing 
effectively performs a homodyne measurement with weak (at the one-photon level) local 
oscillator field. Let us make the substitution ah — > ah + e, where e is the amplitude of the 
mixed drive light, taken for simplicity to be classical, real, and constant. The (unnormalized) 
intensity correlation function is now 

G ( 2 \t,t + r) = G h 2 \t,t + r) 

+e 2 {G h 1 \t,t) + G h 1 \t + r,t + r) 

+2Re[G h 1 \t : t + r) + G^(t,t + r)}} 

+e 4 , (8) 

where G^ft, t + r) = (a k (t)ah(t + r)) is the first-order correlation function of the H mode, 
and G h a \t,t + r) = (a h (t)a h (t + r)) is the if-mode anomalous correlation. Third-order 
correlations vanish for weak drive because scattering a second photon [blue wavy lines in 
Fig. [U^a)] leaves the atom in a manifold orthogonal to that reached after scattering one 
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50 100 150 200 

FIG. 4: Evolution of spontaneously created coherences for one fixed atom; off-diagonal matrix 
elements |(|<7_|_i|p(i + r)|<?_i)| (solid line) and \{g±2\p(t + t)|<7q))| (dashed line) are plotted as a 
function of the time after a first if-photon is detected. The parameters are: 5/7 = 0.25, re/7 = 0.5, 
£/j = 0.3, and a Larmor frequency of 2.27/3. The number of F-photons inside the cavity is 
approximately 0.3. 

photon [red wavy lines in Fig. DJa)], i.e., one- and two-photon states entangle with orthog- 
onal atomic states. The anomalous correlation survives because the manifold reached by 
scattering two photons contains the initial state — \go) in Fig. (Ha). 

With increasing e, the frequency of the quantum beat changes from twice the Larmor 
frequency to the Larmor frequency as the third term on the righ-hand side of Eq. (JHJ) comes 
to dominate the first. Figure [5] illustrates the transition. 

B. Many atoms 

Our experiment is performed with a cold atomic beam and therefore the many-atom case 
must be considered. We generalize the non-Hermitian Hamiltonian of Eq. ([3]) by writing the 
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Fractional drive 



FIG. 5: Size of the beat at the Larmor frequency (solid line) relative to that at twice the Larmor 
frequency (dashed line), as a function of the fraction of drive light mixed with the iJ-mode. The 
parameters are: g/j = 0.25, k/j = 0.5, £ /'y = 0.3, and a Larmor frequency of 2.27/3. 



free, interaction, and non-unitary loss terms for iV arbitrarily located atoms: 
Hq = %u v a\a v + hujha\ a h 

N I 3 4 

+ 2 Z) fit*>tn\9i) 3 (9i\ j + Z fru) ei \ei) j fa 

j=l \i=-3 t=— 4 



N 

Hi = E HHK + 4(Sj + + 4J] + h-c. , (9b) 

H L = —ihK(ala v + a[a h ) 
N «t 

-j2ihi(z$zi + + sitsij, (9c) 

where the dipole coupling constants, gj, j = 1, . . . , N, vary with the location of the atoms 
within the cavity mode function. 

We aim to write the measured correlation function as a sum of terms applying to emission 
pathways for different atoms and their interference. To this end, we formally integrate the 
Heisenberg equation of motion for the if-mode annihilation operator, including its coupling 
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(9a) 



to the reservoir. This yields (see, e.g., Ref. page 206) 

+ EW z{{t-t')e- Kt dt' + v.i., 

j= i J o 



N 

a h {t) = a h (0)e- Ki 



(10) 

E« fc (t) = E*_(t)+E* + (t), (11) 

where v.f. indicates the presence of a vacuum field (reservoir) noise operator, which may 
be dropped from the calculation of a normal- and time-ordered average. We also drop the 
first term on the right-hand side of Eq. fflUj) . as we are interested in times much longer than 
the cavity decay time. Adding the mixed drive amplitude e, as above Eq. (JSJ), the if- mode 
cavity output is treated with the substitution 

N 

a h (t)^J2 A j(t) + £, (12) 
A s {t) = g, I Y? h {t - t')e- Kt 'dt'. (13) 

J 

We further assume that (i) the probability that an emitted if-photon be re-absorbed before 
leaving the cavity is negligible and (ii) no additional if -photon is emitted in between the 
detection of photons at t and t + r. The assumptions are justified, respectively, for moderate- 
to-weak dipole coupling and weak drive. They allow us to treat the atoms as independent 
and write the intensity correlation function as a generalization of Eq. (JHJ): 

N / N 

G®(t,t + T) = £ Gf(M + r)+ E G5 a) (M + r)(Gl a) (M + r) 

3=1 V k?j=l 
N N 

3=1 kjtj=l 
N 

+e 2 E {G? ] {t, t) + Gf (t + r,t + r)+ 2Re (t, t + r) + Gf {t, t + r)]} ^) 

3=1 

where we introduce individual atom correlation functions: 

Gf \t,t + r) = (A](t)A](t + T)Aj(t + r)Aj(t)), (15a) 

Gf\t,t + r) = (A](t)Aj(t + r)), Gf\t,t + r) = (A](t)A](t + r)>. (15b) 

In order to help with the interpretation of this expression, let us assume stationarity 
(correlation functions independent of t) and identical atoms. Equation (TH|) then reduces to 

GW(t) = NG%\t) + N(N-1)\g£\t)\ 2 
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+iV(iV-l)[/i + |Gi 1) (r)| 2 ] 
+2e 2 N{I A + Re[G? (r) + G%\r)]} 

+e\ (16) 
where I A = G A \o) and the subscript denotes any atom. The result is similar to Eq. (4) in 



28j |. The only change is the anomalous correlation G A (t), which is retained here because 
we consider scattering into a single cavity mode, not free-space scattering as in j^sj; in the 
latter case, the random phases accompanying propagation from randomly located atoms to 
the detector cause the anomalous correlation to vanish. It is also important to recall that 
(r) is only nonzero because scattering two iJ-photons places the atom in a ground-state 
manifold that is not orthogonal to the initial state. G^ a '{r) is therefore intimately related 
to the erasure of which-path information by the scattering of a second photon. Of course 



it approaches zero as r — > oo, due to dephasing induced by spontaneous emission [26J]. The 
time scale for this is long compared with the transit time through the cavity if the drive is 
weak. 

Leaving aside the homodyne terms, three correlation functions contribute to the quantum 
beat in Eq. (Tl6l) . There is first and foremost G^\r). It records the beat due to the 
interference of indistinguishable pathways for scattering two photons by one atom; there 
are N such one-atom terms. There are then N(N — 1) copies of the first-order correlation 
function IC^" 1 (r) | 2 . These arise from the sum over cross-terms G^(t)[G^(t)]*, k ^ j, in 
Eq. (|T4l) . They also exhibit a beat at twice the Larmor frequency. It records the interference 
of indistinguishable pathways for scattering a first if-photon from atom j (k) and a second 
from atom k (j) — i.e., the interference of scattering events with reversed time-order, when 
both orders leave the same two atoms in the same final state. Finally, there are N(N — 1) 
copies of (r) | 2 . These anomalous correlation functions record the many-atom extension 

(2) 

of the G A (r) quantum beat. They are present because when one atom scatters two photons, 
our cavity setup is unable to tell from which of the N atoms the two photons come. 

C. Atomic motion 

Our experiment is performed with a slow atomic beam. The atoms move through the 
cavity mode function and the dipole coupling coefficients, gj(t), j — 1, . . . , N, are randomly 
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determined functions of time. If we neglect the effect of photon scattering on the atomic 
motion, gj(t), for a particular atom, is defined by a (constant) velocity Vj, a time, tj, at 
which the atom crosses the plane perpendicular to the beam containing the cavity axis, 
and a position rj on that plane; the velocity may be further specified by a speed Vj, and 
polar and azimuthal angles, 9j and (f)j, defined with respect to the cavity axis. The quantity 
relevant for our experiment is the time average 



The sums in Eq. (TH1) may be considered to range over every atom that enters the cavity 
during the course of the experiment. The atoms are uniformly distributed in a beam of 
rectangular cross section, height d and extension I along the cavity axis, at flux density F 
(number of atoms per second through unit area). Correlation functions for atom j depend 
on Vj, Vj, and the time difference t — tj. They fall to zero for \t — tj\ much larger than the 
mean transit time across the cavity, r = w /(v), where Wq is the mode function waist and 
(v) is the mean speed of an atom along the axis of the atomic beam. 

We consider the mean intensity as an illustration of the way to proceed and then pass 
directly to the result for G^ 2 \t). For the time integral on the right-hand side of Eq. (TTBl) . 
we may write (with t' — t — tj) 



where it is clear that the number of nonzero contributions arising from the sum is just the 
mean number of atoms, FTid, crossing the plane containing the cavity axis during time T, 
while the angle bracket denotes an ensemble average, over rj and Vj, for atoms distributed 
within the beam cross-section; the ensemble average can be taken numerically. We then 
define a mean intensity per (effective) atom normalized to the profile of the mode function, 




(17) 



of Eq. (fl4|) ; we need also, for normalization, the mean intensity 




(18) 




(19) 




(20) 



and arrive at 



/ = N cS I A + e 2 , 



(21) 
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where N e s = pttwq£/4 is the effective number of atoms (see Ref. 29(, for example), with 
p = Fj (v) the atomic density. 

The time average of each sum in Eq. (j!4p is treated in a similar way. This yields a 
straightforward generalization of Eq. f lT5|) : 

G^(r) = N eS G%\r) + N? s \G%\t)\> 
+iV e 2 ff [/i+|G«(r)| 2 ] 
+2e 2 N cS {I A + Re[G^(r) + G%\t)]} 
+e\ (22) 

with (f = 1,2, a) 

G ^=^vSiy^' t,+T)i ")- <23> 

We compare experimental results with a numerical evaluation of Eq. (1221) in Section [III CI In 
the next subsection we introduce simplifications that lead to a closed expression for G^ 2 \t). 



D. Bad-cavity and adiabatic limit 

We make two simplifying assumptions and focus on the weak- field limit. First, we assume 
that the cavity decay rate k is sufficiently large compared with all other rates that the dipole 



operator may be taken outside the integral in Eq. (fT3|) (bad-cavity limit 27|) ; allowing us to 
write Aj(t) = gj(t)Y? h (t) / k. We assume also that the atomic motion is slow and atomic states 
follow the changing coupling constants, gj(t), j — 1, . . . , N, adiabatically. Then noting that 
the dipole operator is expected to carry an atomic excitation proportional to gj(t) in the 
weak-field limit, which we verify numerically, we may factor out the (^(independence and 
write each individual atom correlation function in the form 



G^lv^tjitjt + r) 



G<2x(t), (24) 



K 

with z e = 2,4,2 for f = 1,2, a, where G£ x ( r) is the correlation function, in the long-time 
limit, for a fixed atom at maximum coupling strength g max . The scaled coupling functions, 
gj(t), j = 1, . . . , N, follow from the cavity mode function: 

gj(t) = cos[A;^(t)]e-^ (t)2+ * (i)2]/u; o , (25) 
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where k = 2ti/\ is the wavenumber, and Xj(t), yj(t), Zj(t) define the trajectory of atom 
j. We assign the z-axis parallel to the cavity axis and the x-axis as the axis of the atomic 
beam; at time t' = t — tj = 0, we identify rj = (yj(tj), Zj(tj)), and Vj = (xj(t),yj(t), Zj(t)). 

With these simplifications, the time integration and average of Eq. f[2"3"j) are defined en- 
tirely by the prescribed coupling functions (jj(t), j = 1, . . . , N, and the velocity distribution 
of the atomic beam. In some simple cases analytical results can be obtained. The first 
is when the atoms move along parallel trajectories perpendicular to the cavity and at a 
common speed v. Equations ( 123]) and ( )24|) then yield 

G { 2\r) = ||^f e^^G^Ur), (26a) 

and 

G (l,a) {r) = 3fl* -^fca) (r) (2 g b) 
8 K z 



and setting e = 0, for simplicity, and introducing the normalized correlation functions 

a W( T ) - G(2) ^ (2) ( T ) - G rnL( r ) (2? N 



max 



and = 1, a) 



9«L(r) = < 27b > 



we arrive at a correlation function with Gaussian transit time decay: 



35 



(28) 



Note how the dominant term is 5' m ! ax( r ) (interfering pathways for two photons emitted by 
one atom) at small values of N e s, less than one, whereas for iV e ff 3> 1, l^mLl 1 ")! 2 (interfering 
time orders for two photons emitted by different atoms) and |fi , m I ax( r )| 2 (interfering pathways 
for two photons emitted by one atom or by another) are the dominant terms. 

To illustrate the quantum beat we recover the case of fixed atoms (averaged over locations) 
by taking tq — > oo. Two examples are shown in Fig. E{a): the first with all three terms 
making equal contributions to the beat, and the second where |5 , max( r )| 2 an d |fi , max( T )| 2 
dominate. Figure EJ^b) illustrates the behavior of each term separately for the former case. 
All terms oscillate with the same frequency — twice the ground-state Larmor frequency — in 
the weak-drive limit. There is a phase difference between the oscillations of 9^1 x (t) and 
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FIG. 6: (a) Quantum beat plotted from Eq. ([28]) in the fixed atom limit (to — > oo): iV e ff = 35/36 
(solid) and N c g S> 1 (dashed), (b) Intensity J(r) of the three terms that contribute to make up the 
oscillatory part of the solid curve in (a): gmax(T) (solid), |<?mix(T)| 2 (dashed), and |<?max(T)| 2 (dot- 
dashed). The parameters are: <? m ax/7 = 0.005, k/j = 0.5, £/j = 0.005, and a Larmor frequency 
of 7 /3 



l#max( r )| 2 (two-photon amplitudes) and those of |5 , max( r )| 2 (one-photon amplitudes). This 
arises from different gyromagnetic ratios in the ground and excited states. 

A more realistic modeling of our experiment takes the speed distribution, D(v), of the 
atoms to correspond to a thermal effusive source: D(v) = 2a~ A v 3 e _1,2 / Q2 , where a 2 = 
2kBT/m 30]. The additional average over speed yields 



9 i2) (r) = 1+1 + ^ (bSL(r)| a + b£L(r)| 2 ) 
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+ !«H3~ &{T) ' (29) 

where tq = wq/ (v) is defined with (v) = J^ksT/itm — the mean speed of an atom in the 
source. 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 

x (\is) 

FIG. 7: Effect of a mean tilt of the atomic beam relative to the cavity axis. The envelope of 
(/( 2 )(t) — 1 is plotted for no tilt of the beam (solid) and a mean tilt (6) = 1.3deg (dashed), with 
(v) = 22m/s, Av = 2m/s, and AO = 0.025 deg. 

In practice the atoms do not move perfectly perpendicular to the cavity axis. We do not 
have an analytical expressions for this most general case. Figure [7] shows the numerically 
calculated envelope of the correlation function for the case N c g 3> 1, assuming a triangular 
distribution for the polar angle relative to the cavity axis 9 and a Maxwell-Boltzmann speed 
distribution. The figure shows how the standing-wave structure of the cavity mode function 
becomes imprinted on the envelope of the quantum beats. The local minimum of the dotted 
curve corresponds to the delay time when a majority of atoms pass from an anti-node to a 
node between the detection of the first and second photons. The spread in angle and speed 
explains why this structure does not recur at longer delays. 
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III. EXPERIMENT AND RESULTS 



A. Apparatus 

We perform measurements using a slightly modified version of the apparatus described 
in Ref. {31]. A sketch of the experimental setup appears in Fig. QJb). We probe a small 
ensemble of 85 Rb atoms coupled to a Fabry- Perot resonator in vacuum. The 2.2 mm cavity 
has a 56 /im mode waist and a finesse of 11,000, with losses shared approximately equally 
between mirror transmission and scattering/absorption. The decay rates of the field and 
atomic dipole — (k,j/2)/27T = (2.8,3.0) x 10 6 s~ 1 — are approximately matched, and twice as 
large as the dipole coupling strength, g max /27i = 1.5 MHz, on the D 2 -line F = 3, m = to 
F' — 4, m' = transition; this places our system in the intermediate coupling regime of 
cavity QED — single-atom cooperativity C\ = g 2 /jk = 0.12 and saturation photon number 
no = 7 2 /3g 2 = 5.3 — with only a small probability of reabsorption after a photon is emitted 
into the cavity mode. 

A crossed-polarizer configuration separates the weak if -mode fluorescence from the much 
stronger V-polarized drive, necessitating careful selection and alignment of polarization el- 
ements. We drive the cavity with a laser sideband generated by a polarization-maintaining 
single- mode 780 nm fiber modulator (EO Space) operated at 230 MHz. Before entering the 
mode-matching lens and vacuum chamber, the drive is linearly polarized — extinction ra- 
tio less than 5 x 10~ 5 — after passing through a Glan-Thompson polarizer and zero-order 
half- wave plate (HWP). A second zero-order HWP placed after the cavity aligns the po- 
larization to a calcite Wollaston prism for separation of the H- and V-mode light at the 
output. The extinction ratio after this splitter is limited by birefringence in the cavity mir- 
rors, vacuum chamber windows and lenses. Its exact value is a function of drive intensity, 
likely due to thermoelastic stress-induced birefringence in the components, but is generally 
of order 5 x 10~ 4 . The splitting of H- and V-mode resonance frequencies due to cavity 
mirror birefringence is less than 200 kHz. 

The separated beams go to two avalanche photodiodes (APDs, Perkin-Elmer) for photon 
counting, except when measuring autocorrelation functions, in which case the V-mode beam 
is blocked and the if -mode beam is split equally between the two detectors by means of 
a separate HWP and polarizing beam splitter. The TTL output channel of each detector 
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is electronically split between a counter unit, for measuring rates, and a PC time-stamp 
card (Becker and Hickl DPC-230) for recording detection times with 164ps resolution. A 
typical time series measurement takes approximately 300 s. The output of the time-stamp 
card is written to a text file and parsed by a C++ program to calculate the cross-correlation 
between detection events (the autocorrelation of the if -mode). The correlated events are 
recorded in a histogram, bin width 16.4 ns, extending out to ±16.4 /is. Division by the 
average (uncorrelated) bin count yields g^ 2 \r). The power spectrum is calculated from a 
discrete fast fourier transform (FFT) of this function. 

The 85 Rb atoms are extracted continuously from a magneto-optical trap (MOT) operating 
as a low-velocity intense source (LVIS) 32j in a chamber directly above the cavity. The 
atoms have a mean speed of ~ 22m/s, which yields an interaction time of a few microseconds. 
Atoms leaving the MOT are primarily pumped to the m = 3 ground state, with quantization 
axis provided by the residual vertical magnetic field from the MOT "anti-Helmholtz" coil pair 
(~ 7G at the location of the cavity). In order to change the magnetic field in the cavity, a 
third coil is added directly below the "anti-Helmholtz" pair, with separate currents applied 
to each of the three coils. We are thus able to vary the vertical magnetic field in the cavity 
between ±12 G while maintaining the required gradient for the MOT. An additional pair 
of coils is oriented with axis parallel to the cavity axis to cancel any residual field in that 
direction. 

Before atoms enter the cavity mode they are optically pumped to the m = ground state 
using a beam resonant with the F = 3 to F' = 3 transition and polarization parallel to the 
vertical magnetic field. The optical pumping beam is combined with light from the MOT re- 
pumper laser in a 50/50 polarization- maintaining single-mode fiber splitter, and collimated 
to a waist diameter of 0.5 mm in order to pass between the top of the cavity mirrors and 
the upper edge of the vacuum window. Due to strong scattering from multiple reflections 
into the APDs, we are unable to use the beam in a retro-reflected configuration. It therefore 
imparts a net momentum kick to the atoms. The intensity of the optical pumping beam is 
chosen optimally as a compromise between moving most atoms to m = 0, while not ejecting 
too many from the beam, and scattering too strongly into the APDs. 

Figure [8] shows a typical measurement sequence used for optimizing the optical pumping 
configuration. Frame (a) shows the absorption fraction (output intensity over input inten- 
sity, denoted X/Y) as measured from the V-mo&e count rates with the cavity and drive 
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FIG. 8: (Color online) (a) Example of the T^-mode absorption fraction as a function of common 
laser and cavity detuning from the central m = to m! = atomic transition; for a magnetic 
field of ~ 5G and approximately one photon on average in the cavity. The optical pumping beam 
(red/open circles "on", blue/filled dots "off") shifts the population toward the center. Solid lines 
are least-square fits of the data to a Lorentzian line shape, (b) if-mode count rates measured 
concurrently with the data in (a). 



frequencies simultaneously swept across the atomic resonance. The effect of the optical 
pumping is to shift and narrow the absorption peak, ideally to yield a symmetric lineshape 
centered around the resonance frequency of the F = 3, m = to F' = 4, ml = transition 
(0MHz in the plot). The departure of the center frequency from zero results from a combi- 
nation of incomplete optical pumping and a small drift in the frequency setpoint of the Rb 
saturated absorption spectroscopy reference used for the laser. In frame (b) count rates for 
the if- mode light show a similar effect. We use the center frequency obtained from these 
scans as the reference (zero-detuning) point for our measurements. 

The extensive Optical Bistability literature is useful for understanding and interpreting 
the effects of absorption and detuning in our experiment. When making these connections 



(see for example 



331 ]) it should be noted that we operate in the low intensity limit, and 



simultaneously scan the laser and cavity in order to address the atoms directly. 



B. Results 

Frame (a) of Fig. [9] displays a measured correlation function for a magnetic field of 5G; 
its calculated power spectrum is displayed in frame (b). The main peak near 4.8 MHz 
corresponds to the quantum beat resonance. A smaller peak at half this frequency is also 
present, though largely obscured by noise. It is the result of homodyne interference with 
drive light mixed in by cavity birefringence [see the paragraph surrounding Eq. (jSJ)]. The 
small sidebands on the main peak correspond to a slight modulation of the beat envelope. 
The modulation is visible in frame (a) and results from the small (1 — 2 degree) deviation 
of the atomic beam from normal incidence with the axis of the cavity, which introduces 
sinusoidally varying coupling strengths, gj(t), and amplitude modulation of the spontaneous 
emission rate (see Fig. [7j). 

Figure fTUl illustrates the changing frequency of the quantum beat with increasing magnetic 
field, where the expected linear dependence is observed. We note that the beat frequencies 
also depend on the intensity of the drive through an anomalous light-shift, which we report 



elsewhere 26|. Those presented here are extracted as the zero-intensity (i.e. unshifted) limit 
of the measured frequencies for each magnetic field. Our zeroing of the magnetic field in all 
three directions is not better than 10 mG. This error in independent calibration makes the 
small offset in the figure consistent with zero. 
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FIG. 9: (a) Measured intensity correlation function g^ 2 \r) and (b) its FFT power spectrum; for a 
5G magnetic field, iV e g = 2.9, and approximately 6.5 photons in the V mode with no atoms present. 
The peak in the spectrum is located at ~ 4.8 MHz, twice the ground state Larmor frequency for 
85 Rb in a 5G field. 
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FIG. 10: Measured linear dependence of the frequency of the quantum beat on magnetic field. 



Figure [IT] illustrates the change in the observed beat when the polarization presented 
to the detector is not taken orthogonal to that of the drive but allowed to rotate by a few 
degrees. The rotation is controlled by changing the angle of the HWP placed between the 
cavity and the PBS (Fig. [JJd). This mixes a small amount of drive light with the scattered 
light. With increasing mixed- in fraction the beat is eventually dominated by a homodyne 
term, which arises from the correlation of a photon scattered into the H mode with a photon 
from the drive [see Fig. [5] and terms proportional to e 2 in Eqs. (JHJ), (tHJ), ( fl6l) . and (122]) ]: 
thus, as in the two-atom case, interfering time orders also yield a quantum beat. This 
beat oscillates at half the frequency and allows the correlation function to dip below one. 
Generally, some drive light is coupled into the H mode through a small birefringence of 
the cavity mirrors. Frames (c) and (d) of Fig. [11] are recorded at the HWP angle that 
gives maximum visibility of the half- frequency beat. Frames (b) and (d) show the dramatic 
increase in visibility gained when the driving laser is slightly detuned from resonance. This 
is related to a decreased decoherence rate from quantum jumps (see Ref. 26|). 

In Fig. [12] we show how g^ 2 '(r) evolves as the number of effective atoms increases from 
less than one to nearly three. For the fewest atoms [frame (a)], fluctuations in the number 
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FIG. 11: Evolution of the measured </ 2 ^(t) with homodyne mixing and detuning of the drive: (a) 
on resonance, maximum extinction; (b) —6 MHz detuning, maximum extinction; (c) on resonance, 
HWP rotated by 2.8 degrees; (d) — 6MHz detuning, HWP rotated by 1.2 degrees. Data taken for 
a 5G magnetic field, iV e g = 2.9, and approximately 2.0 photons in the V mode with no atoms 
present. 
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interacting with the cavity show up as a broad Gaussian background peak, reflecting the 
increase in the scattering rate when an atom is present (similar to Ref. 3l|.) The beats 
sit on top of this background, with small visibility, a consequence of optical pumping and 
spontaneous emission (to modes other than the cavity mode). The correlation function is 
dominated by the contribution from g^l x (r) in Eqs. (I2"gj) and (|2"9"|) . As the density of the 
atomic beam grows [frames (b) and (c)] the background peak disappears as contributions 
from multiple atoms become more prominent; contributions from \g^l K (T)\ 2 , |<7max( r )| 2 ' an d 
fl , max( r ) contribute with more-or-less weight in frame (c). 

Finally, Fig. [13] shows the change in the spectrum when the optical pumping beam is 
added. The effect on the beat frequency is minimal, revealing the robustness of the quantum 
interference to the initial distribution of the atoms amongst the ground-state Zeeman levels. 
The biggest change is in the low-frequency components of the FFT, which correspond to 
atomic motion through the standing-wave mode. This is caused by the momentum kick 
imparted by the optical pumping beam. 



C. Comparison of theory and experiment 

Outside the bad-cavity and adiabatic limit, the normalized correlation function is given 

by 

g^\r) = G^\r)/l\ (30) 

with / and G (2) (r) defined in Eqs. (ED and ([22]). For comparison with experiment, we 
introduce a global scale parameter, s, writing 

9 (2 \r) = l + s^\r)\ 2 + \g^(r)\ 2 

+ Areff f +e2 2Re[^ 1) (r)+^ ) (r)]| , (31) 



where (£ = 1, 2, a) 



iVeffi + e z 



This expression is evaluated numerically by calculating the one-atom correlation functions 
from the model of the atomic beam developed in Section III CI The speeds Vj are selected 
from a Gaussian distribution, and the angles 0j and <pj from a triangular distribution. 
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FIG. 12: Evolution of the measured g^ 2 \r) with increasing atomic beam density. Data taken for a 
5G magnetic field, iV e g = 0.3,0,9,2.9, (a,b,c) and approximately 3.1 photons in the V mode with 
no atoms present. 



Figure ITffi a) and (b) show a fit to the experimental data and its FFT power spectrum. 
Experimental error bars are computed as the square root of the number of photon counts 
in each bin. In order to improve the fit, we adjust the mean speed and angles according to 
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FIG. 13: (Color online) Sample power spectrum [FFT of the measured g( 2 \r)] with (red, dashed) 
and without (blue, solid) optical pumping of the atoms prior to entering the cavity. 



the following considerations. The experimental data shows a low frequency modulation at 
short delays, the signature of a small mean inclination of the atomic trajectories away from 
normal to the cavity axis (see Fig. [7]). Comparing the FFT of the measured autocorrelation 
function with simulations at different mean beam angles (6), with (0) = 0, we find that 
(9) = 1.4 degrees optimizes the fit. For small angles like this, the main parameter affecting 
the width of the peak at twice the Larmor frequency, around 5 MHz in the figure, is the mean 
atomic speed. A value of (v) = 17m/s optimizes the fit, consistent with values expected 
from an LVIS 32]. The amplitude of the peak around 2.5 MHz is determined by the value 
of e, which is adjusted as another free parameter. 

Frame (a) of Fig. [TH compares theory and experiment in the time domain. The reduced 
chi-square value for this fit is 1.3. The major differences between the calculation and the 
experiment comes for delays close to zero; here the measured correlation function is substan- 
tially larger in value compared to the calculated one. This is at least partially due to the 
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presence of uncorrelated background light scattered into the detection path (primarily from 
the MOT cooling beams). Frame (b) of Fig. [T41 shows the fit in the frequency domain. The 
calculated spectrum in fact goes to zero at high frequencies (larger than 7 MHz), while the 
experimental one remains flat due to the presence of background light in the detectors (up 
to about 200 MHz). With a flat frequency background added to the calculation to account 
for the residual background and shot noise (as shown in the figure), the reduced chi-square 
value of the fit is 0.99. Frame (c) shows a similar correlation function after addition of a 
small amount of coherent drive in order to enhance the homodyne signal. The simulation 
accurately captures the pronounced change in frequency and shape of the signal, with a 
reduced chi-square value of 1.6. 



IV. SUMMARY AND CONCLUSIONS 



We have studied theoretically and experimentally how the ground-state quantum beats 
reported in 12j depend on different parameters. The fundamental beat frequency occurs at 
twice the Larmor frequency and is found to increase linearly with magnetic field as expected. 
Mixing of the driving and scattered fields produces a beat at the Larmor frequency itself. 
Increasing the number of atoms brings the minimum of the oscillation to the shot noise level; 
this is because the many pairs of two-atom beats come to dominate the one-atom signal. 
The stochastic evolution of coherence within the atomic level structure shows that the many 
levels of the F = 3 to F' = 4 transition help make the observed quantum beats robust 
against optical pumping. 



Appendix A: Dipole operators 

The one atom dipole operators are given by 



s - = J \\9-3)(e- 3 \ + J^\g- 2 ){e- 2 \ + J^j|#-i)(e_i| + J^\g )(e \ + J^\g+i)(e +1 \ 



+y yb+2><e +2 | + J^\g +3 )(e +3 \ 



(Al) 



- - \/ 28 l#-3>(e- 2 | + J^|0_ 2 )(e_i| + \l^\g-i)(e \ + J^\g )(e +1 \ + J^\g +1 )(e +2 \ 
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FIG. 14: (Color online) Sample comparison of the measured (points) and calculated (solid line) 
correlation function, (a), and its FFT power spectrum, (b), under conditions where atomic spon- 
taneous emission dominates the measured photon counts. A flat background is added to the 
calculated spectrum in order to account for residual noise in the detectors. The atomic beam 
parameters are: iVeff = 3, s = 0.53, (v) = 17m/s, Av = 2m/s, (6) = 1.4 deg, AO = 0.9 deg, 
((/>) = deg, Acj) = 0.7 deg. The mean number of scattered photons in the V mode is 0.63, and 
the birefringence background e 2 is 1.2% of the driven V-mode photon number, (c) Sample mea- 
sured (dashed line) and calculated (points) correlation functions with a small amount of coherent 
drive admixed in order to enhance the homodyne signal. Parameters are: N e g = 0.55, s = 1, 
(v) = 13.5m/s, Av = 2m/s, (9) = 0.97 deg, Ad = 0.9 deg, (4>) = deg, A(f> = 0.7 deg. The mean 
number of scattered photons in the V mode is 1.2, and the mixed coherent field e 2 is 0.5% of the 
driven y-mode photon number, which due to the smaller number of atoms is strong enough to 
substantially change the signal. 



33 



+J -\g +2 )(e +3 \ + \g+a) (e+4 1 , 



(A2) 



/3 /15 / 5 / 3 

-b-2>(e_ 3 | + ]l — \g- 1 )(e_ 2 \ + J — \g )(e-i\ +]l—\g + i)(e \ 

3 l~T~ 

'^l#+2>(e+i| + W— |^+3)(e+ 2 | • 



(A3) 



Acknowledgments 



Work supported by NSF, CONACYT, Mexico, and the Marsden Fund of the Royal 
Society of New Zealand. 



[1] U. Eichmann, J. C. Bergquist, J. J. Bollinger, J. M. Gilligan, W. M. Itano, D. J. Wineland, 

and M. G. Raizen, Phys. Rev. Lett. 70, 2359 (1993). 
[2] W. M. Itano, J. C. Bergquist, J. J. Bollinger, D. J. Wineland, U. Eichmann, and M. G. Raizen, 

Phys. Rev. A 57, 4176 (1998). 
[3] G. S. Agarwal, J. von Zanthier, C. Skornia, and H. Walther, Phys. Rev. A 65, 053826 (2002). 
[4] M. O. Scully and K. Driihl, Phys. Rev. A 25, 2208 (1982). 
[5] E. B. Aleksandrov, Opt. Spectrosc. 17, 522 (1964). 

[6] J. N. Dodd, W. J. Sandle, and D. Zissermann, Proc. Phys. Soc. 92, 497 (1967). 
[7] G. Breit, Rev. Mod. Phys. 5, 91 (1933). 

[8] Z. Ficek and S. Swain, Quantum Interference and Coherence, Springer Series in Optical Sci- 
ences (Springer- Verlag, Berlin, 2005). 

[9] W. W. Chow, M. O. Scully, and J. O. Stoner Jr, Phys. Rev. A 11, 1380 (1975). 
[10] R. M. Herman, H. Grotch, R. Kornblith, and J. H. Eberly, Phys. Rev. A 11, 1389 (1975). 
[11] A. G. Zajonc, Phys. Lett. A 96, 61 (1983). 

[12] D. G. Norris, L. A. Orozco, P. Barberis-Blostein, and H. J. Carmichael, Phys. Rev. Lett. 105, 
123602 (2010). 

[13] A. T. Forrester, R. A. Gudmundsen, and P. O. Johnson, Phys. Rev. 99, 1691 (1955). 
[14] S. Gerber, D. Rotter, L. Slodicka, J. Eschner, H. J. Carmichael, and R. Blatt, Phys. Rev. 
Lett. 102, 183601 (2009). 

34 



[15] E. B. Alexandrov, M. P. Chaika, and G. I. Khvostenko, Interference of Atomic States 

(Springer- Verlag, Berlin, 1993). 
[16] A. Aspect, J. Dalibard, P. Grangier, and G. Roger, Opt. Commun. 49, 429 (1984). 
[17] R. T. Willis, F. E. Becerra, L. A. Orozco, and S. L. Rolston, Phys. Rev. A 82, 053842 (2010). 
[18] M. Schubert, I. Siemers, R. Blatt, W. Neuhauser, and P. E. Toschek, Phys. Rev. A 52, 2994 

(1995). 

[19] A. Aspect, E. Arimondo, R. Kaiser, N. Vansteenkiste, and C. Cohen-Tannoudji, J. Opt. Soc. 

Am. B 6, 2112 (1989). 
[20] J. Javanainen, Europhys. Lett. 17, 407 (1992). 
[21] A. K. Patnaik and G. S. Agarwal, Phys. Rev. A 59, 3015 (1999). 
[22] M. Kiffner, J. Evers, and C. H. Keitel, Phys. Rev. Lett. 96, 100403 (2006). 
[23] M. V. G. Dutt, J. Cheng, B. Li, X. Xu, X. Li, P. R. Berman, D. G. Steel, A. S. Bracker, 

D. Gammon, S. E. Economou, et al., Phys. Rev. Lett. 94, 227403 (2005). 
[24] S. E. Economou, R. B. Liu, L. J. Sham, and D. G. Steel, Phys. Rev. B 71, 195327 (2005). 
[25] H. J. Carmichael, An Open Systems Approach to Quantum Optics, Lecture Notes in Physics, 

vol. 18 (Springer- Verlag, Berlin, 1993). 
[26] D. G. Norris, A. D. Cimmarusti, L. A. Orozco, P. Barberis-Blostein, and H. J. Carmichael, 

Phys. Rev. A 85, 021804 (2012). 
[27] H. J. Carmichael, Statistical Methods in Quantum Optics 2: Non-Classical Fields (Springer- 

Verlag, Berlin, 2008). 

[28] H. J. Carmichael, P. Drummond, P. Meystre, and D. F. Walls, J. Phys. A: Math. Gen. 11, 
L121 (1978). 

[29] H. J. Carmichael and B. C. Sanders, Phys. Rev. A 60, 2497 (1999). 

[30] N. Ramsey, Molecular beams, Oxford classic texts in the physical sciences (Clarendon Press, 
Oxford, 1985). 

[31] D. G. Norris, E. J. Cahoon, and L. A. Orozco, Phys. Rev. A 80, 043830 (2009). 

[32] Z. T. Lu, K. L. Corwin, M. J. Renn, M. H. Anderson, E. A. Cornell, and C. E. Wieman, Phys. 

Rev. Lett. 77, 3331 (1996). 
[33] L. A. Lugiato, in Progress in Optics, edited by E. Wolf (North-Holland, Amsterdam, 1984), 

vol. XXI, pp. 69-216. 



35 



